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In the inflationary scenario driven by a single self-interacting scalar field, the 
dominant mode of primordial fluctuations is predicted to be Gaussian to a very good 
approximation. However, it has long been appreciated that non-trivial information 
should be encoded in the connected three-point function and higher connected 
correlation functions p U 0 II U- These correlation functions are expected to 
leave signatures in the statistical properties of the cosmic microwave background 
(CMB) temperature anisotropies, and may already have been detected [HI ZJ ;S|. In 
principle, such information yields crucial insight into the nature of the underlying scalar 
field Lagrangian during the inflationary epoch, and provides a sensitive discriminant 
jMnmn between the large number of competing inflationary models nansmunanni 

In the standard single field inflationary scenario, the scalar field action is generally 
taken to be of the form 

s = ~f d4 x + , ( 1 ) 

where the kinetic term is canonically normalized. Choosing a suitable inflationary 
model then corresponds to engineering an appropriate form for the potential, V(<j)). 
Unfortunately, identifying an acceptable form for the inflationary potential has proved 
to be a difficult task DU- 

The scalar field potential is not necessarily the only degree of freedom in inflationary 
model building. Indeed, in models descending from a supergravity or superstring 
compactification, where the inflaton might be identified with a light moduli field, it is 
generally expected that corrections to the kinetic term of the scalar field action (HD will 
arise [TSJ. Moreover, even if the description of microphysical degrees of freedom given 
by the action 0 is appropriate at the classical level, one would generally expect loops 
in the quantum theory to generate operators in the Lagrangian that are proportional 
to higher derivatives (<90) 2 , (<90) 4 , and so on. Such interactions would be suppressed 
by powers of the renormalization scale M, and if this scale is large, M ~ Mp, where 
Mp is the Planck mass, the contribution of these operators would be negligible at the 
energy scale of inflation. On the other hand, if M is closer to the unification scale of 
some Grand Unified Theory, ultraviolet corrections of this type might be significant and 
of considerable relevance in the very early universe using. 

Such non-minimal choices of the scalar field action can be written in the form 

S = |d 4 ^P(^,%^,...), (2) 

where P is an arbitrary Lorentz-invariant polynomial of 0 and its derivatives. This form 
of the action includes the standard choice (JTJ) as a special case. Non-trivial choices of 
kinetic terms of the sort described by © have been considered previously by a number 

of authors PU121 HQ 121 mi 121 ESI • 

The presence of higher-derivative operators in P allows for qualitatively new effects. 
For example, in the case where P is a function only of first derivatives of the field (and 
is independent of the value of the field itself), it is possible for the inflaton to ‘condense’ 
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at a turning point in P. In this scenario, the field momentum acquires a non-zero 
vacuum expectation value, (</>) 7 ^ 0 [23 , that is able to drive a phase of de Sitter 

(exponential) expansion. Low-energy approximations to stringy degrees of freedom may 
also be described through an action of the form ©>• Specific examples include the 
tachyon field or a generalized Dirac-Born-Infeld action |2S11123- 

in view of the above possibilities, therefore, a study of more complicated actions 
of the form © is well motivated. In general, one may include terms with any number 
of higher derivatives in P. However, in any effective theory it is to be expected that 
terms containing higher derivatives will be suppressed by powers of the ultraviolet cut¬ 
off scale. In theories coupled to Einstein gravity, this is most naturally the Planck 
scale, so if inflation occurs at energies E <C Mp, the contribution from such operators 
eventually becomes small. This implies that the effect of arbitrarily higher derivatives 
can be neglected. We therefore consider the class of theories where P contains at most 
first derivatives in the scalar field 

P = P(X, 0), X = -<f 6 V a 0V 6 0. (3) 

Having restricted the system to first derivatives in this way, the requirement that these 
derivatives enter via A" is fixed by Lorentz invariance, so this really is the most general 
form of the Lagrangian. 

Despite the novel character of ©, and the widely differing physics that enters 
into the various models leading to such an effective Lagrangian, the predictions for 
standard observables such as the scalar spectral index are essentially degenerate with 
the standard scenario o to leading order in the slow-roll parameter e ©EES!, where 
e = - H/H 2 < 1 . This implies that further observational information is required 
in order to discriminate between the alternative scenarios. One source of additional 
observational insight is provided by the three-point or higher connected correlation 
functions, as measured in principle through CMB fluctuations. Since the largest 
contribution is expected to arise from the three-point function, considerable attention 
has focused recently on the theoretical nature of this correlation and its possible 
observational detection in the CMB 21] • 

In this paper, we calculate the three-point function for a general theory of single¬ 
field inflation whose microphysics is described by the action ©-©, using the slow-roll 
approximation to control the calculation where necessary. The three-point function 
for the canonical action © coupled to gravity was calculated by Maldacena i.22j ( see 
also [3]), and a similar analysis has been performed by Rigopoulos & Shcllard |28j . 
who approximate the quantum fluctuations by a stochastic noise term. Some specific 
examples of theories containing higher-derivative operators have been considered in the 
literature H3 ESI 112 ESI, either coupled to gravity or considered in isolation. We 
perform the computation quite generally, including all details of the minimal coupling 
to Einstein gravity. Since non-gaussianity is a potentially sensitive probe of new or 
unexpected physics, some attention has also been given to non-standard scenarios, such 
as tachyon or brane inflation [29 . The case of scalar field inflation with canonical kinetic 
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terms coupled non-minimally to gravity was recently considered by Koh, Kim & Song 

EDI- 

One of the key features that arises in considering the generalized action 
is that the speed of sound, c s , is in general time-dependent and differs from unity (in 
units where the speed of light c = 1 ), in contrast to the canonical action (fTl) . where 
c s = 1 . As we shall see, this has significant implications for the form of the three-point 
correlation function. Our main result is that the three-point function contains terms 
that have similar ^-dependences to that of standard, single-field models ra, but with a 
different dependence on the slow-roll parameters. Furthermore, new fc-dependent terms 
are present which are entirely absent in the standard case and arise whenever the speed 
of sound differs from unity. In principle, these new features represent a distinctive 
probe of c s 7 ^ 1 in the CMB [[101, since they produce a pattern of angular dependence 
on the CMB sky which ought to be accessible whenever the non-gaussianity predicted 
by M is observable. Our expression for the three-point function also respects the 
long-wavelength gravitational consistency relation 1221313. 

The outline of this paper is as follows. We discuss the homogeneous background 
model in Section [21 establishing our notation for the equations of motion and proceeding 
to discuss the slow-roll approximation for models with generalized kinetic terms. In these 
models, there is an extra requirement, over and above the familiar restriction that held 
derivatives should be less than the expansion rate, which follows from demanding that 
the rate of change of the sound speed should be sufficiently small. In Section |3] we 
couple fluctuations in the scalar held to scalar modes of the metric. This is most simply 
expressed using the ADM decomposition [32] , We solve the ADM constraint equations 
in Section IdTl Given a solution of the constraints it is possible to construct the Gaussian 
action, as shown in Section Id.21 This has been done previously in the literature using 
different techniques 1201122 - 

In Section El we construct the interaction vertex for the coupled huctuations. This 
interaction vertex may be viewed as the generalization of the third-order action presented 
in to include a varying speed of sound, or alternatively, as an extension to third- 
order of the fc-inhation action derived by Garriga & Mukhanov |20j . We calculate the 
three-point function by employing a different technique to that most commonly used 
in the literature to date. This provides an alternative and independent method for 
calculating the three-point function which may sometimes prove more convenient. We 
outline this method in Section EJ and show how it leads to an expression for the tree-level 
(‘semiclassical’) approximation of the three-point function in a particularly transparent 
and direct way. The details of the calculation are presented in Section El We recast 
the result as an expression for the conventional non-linearity parameter /nl and show 
explicitly how the result obeys Maldacena’s consistency condition [22- Finally, we draw 
our conclusions in Sectional 
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We work with an action of the form 

S = ^ J d 4 x [R + 2 P(X, 0 )], X = -g ab X a fX b f, (4) 

where units are chosen such that the reduced Planck mass Mf 2 = 8 itG is set to unity. 
The homogeneous background solution is assumed to be of Friedmann-Robertson- 
Walker (FRW) form with flat spatial slices, 

ds 2 = — df 2 + a 2 (t) dx 2 . (5) 


Given this metric, the gravitational dynamics supply both an evolution equation and a 
constraint for the Hubble parameter, p[ = a/a: 

2 H + 3 H 2 = - P , 

H 2 = i( 2 A :P,x - P), ( 6 ) 

where a comma denotes a partial derivative. The constraint equation is merely the 
Friedmann equation in this model. The equation of motion for the scalar field is given 
by 

X(P,x + 2 XP, xx ) + 2V3(2 XP x - P)^ 2 XP x = X l ' 2 (P^ - 2 XP X(j> ). (7) 


An important consequence of the non-trivial kinetic structure in P is that the naive 
dispersion relation for 0 is modified, and fluctuations in the scalar field do not travel at 
the speed of light. Instead, the sound speed in 0 is given by 
2 _ P,x 

Cs ~ P, x + 2XP xx ' 


2.1. The slow-roll approximation 

For general P the scalar field equation CD cannot be solved analytically. In order 
to proceed, therefore, it is necessary to resort to approximations, where the solution is 
expanded perturbatively in powers of a small parameter. Within the context of standard 
scalar field inflation, this is usually achieved by assuming that the field 0 is rolling slowly 
in comparison to the expansion rate, i.e., 0 2 -C H 2 . More quantitatively, we may define 
H _ XP iX _ e 

s ~ ~1 p ~ h 2 ’ v ~ 7 h' 

with the understanding that |e|, |//| <C 1 for reliable calculations!. In practice we will 
assume that e ~ r/ and express this condition by writing e, rj 0(e). 

f In standard, single-field inflation, e is usually a positive quantity by definition. That need not be 
the case here. There are several inequivalent definitions of g which are used in the literature, of which 
the most common alternatives to our choice are gv = V"/V i<21 an d 1H = 2H"/H JJJl. The former 
definition makes sense only for standard inflation, whereas the latter can be used where non-trivial 
kinetic terms are present. With a canonical choice of kinetic term, one can show that these alternatives 
are related to our g by the rules 

g = -2g H + 2e = - 2g v + 4e. 

Note that these only apply for standard inflation. 


( 10 ) 




Primordial non-gaussianities in single field inflation 


6 


It proves useful to decompose the parameter £ into two new dimensionless ratios, 
£$ and £x, which measure how the expansion rate varies with the kinetic and potential 
parts of 0 , respectively: 


f dH XdH 

£ ~ ~ Ipdx ~ £<p + £x ' 

The scalar held equation of motion 0 may then be written as 


( 11 ) 


X = - 6 Hc 2 s X (l - , (12) 

and this allows us to express t-derivatives in terms of derivatives with respect to A". 

In principle, there is no requirement from a dynamical point of view that £$ and 
£x should both be small, even when |e| <C 1. In standard, single-held inflation, £ 
and g are often referred to as the slow-roll parameters, and the limit |£|, \q\ <C 1 as 
the slow-roll limit. This terminology is not quite appropriate for a general choice of 
P, since |f| -C 1 no longer entails f 2 -C H 2 . For brevity, however, we can refer to £ 
and 77 as flow parameters, since they describe how the theory evolves on the space of 
inflationary models jJS ESj- By an abuse of terminology, we will continue to describe 
the limit |£| <C 1 as ‘slow-roll’ because the content of the approximation is familiar in 
the literature. 

As well as the familiar conditions |£|, \rj\ -C 1, it will also be necessary impose 
bounds on the rate of change of the sound speed due to the generalized kinetic terms in 
0 m EU 231122] • We therefore define the parameters 


u = 



-2A 


P.xx 

T7’ 



(13) 


where s represents a dimensionless measure of the rate of change of the sound speed 
c 2 = 1/(1 — u). These two parameters are related by 


u = 2Hs(l — u ). 


(14) 


It is well-known that the time derivatives of £ and 77 are second-order in the slow-roll 
expansion, in the sense that s, ?) ~ 0(£ 2 ) 34]. This means that we can consistently work 
to first-order in 0(f), while keeping £ and 77 constant. Eq. (fT41) implies that s is related 
to the time derivative of u, so it is sufficient that u = 0(f) in order that s = 0(f 2 ). In 
this case, c 2 departs from unity only by a quantity that is first-order in slow-roll. 

After combining m with the scalar held equation m, we may write down a 
relationship between the parameters £, ?/ and u\ 
f( 2 f - 77 ) - 6£ X 


u = 


f( 2 f - 77 ) - 3f x 

It follows that a necessary condition for u = 0(f) is that £x satisfies 
6£ X = e( 2 f - 77 ) + 0 (f 3 ) 


(15) 


(16) 


and this implies that £x is subdominant with respect to f^. Indirectly, this is a rather 
non-trivial condition on P(A, <f ) and means that no guarantee can be given that a 
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particular P(X,(ft) will necessarily support a phase where |e| -C 1, even in principle. 
In this paper, we do not attempt to ascertain the conditions under which a particular 
P(X, (ft) will admit a slow-roll epoch, but merely provide an expression for the three- 
point function which is valid whenever it does. 

In order to simplify some of the formulae that follow, it will prove useful to introduce 
two new quantities, £ and A, which are combinations of derivatives of P, and defined 
by 


£ = XP )X + 2 X 2 P, 


,XXi 


A = X 2 P xx + -X 3 P xxx . 


These can be written in terms of flow parameters: 


£ = 


H 2 e 


= H 2 e(l-u), 


A = 


2 £ 
3 £ X 


(1 — u)s — u 


(17) 


(18) 


3. The ADM formalism 


Any consistent cosmological calculation of fluctuations in some scalar held (ft which 
dominates the energy density of the universe must account for the universal coupling to 
gravity, since any perturbation in (ft will produce a non-negligible perturbation in the 
energy-momentum tensor. Thus, we need to calculate the action for small fluctuations 
around the homogeneous background solution of (0J), taking into account both the 
perturbations in the scalar held, deft, and the scalar modes of the metric. There is no 
need to include vector perturbations, which die away rapidly with the cosmic expansion 
and are not sourced by inflation. In addition, we omit tensor modes. In principle, tensor 
modes corresponding to gravitational waves are excited by inflation, but gravitational 
waves have not yet been detected and it is anticipated that any non-gaussianity involving 
such modes will be at a lower level than that predicted for the scalar sector [?71 . In 
the near future, observational effort is likely to be directed towards the determination 
of the scalar non-gaussianity, to which we restrict our attention. 

An arbitrary scalar perturbation of the background © can be written in the form 

ds 2 = -(1 + 2<h)dt 2 + 2 a 2 {t)B, l d x'dt + a 2 (t) [(1 - 2T)% + 2 E tij \ dx\ixf (19) 

where a comma denotes a partial derivative with respect to the spatial coordinates 
x l . One could directly calculate the action for the helds <h, B, T and E and work 
with these fluctuations together with fluctuations S(ft in the inflaton. After integrating 
by parts, dropping total derivatives, applying the constraint equations and using the 
background equations of motion, it can be shown that to quadratic order the action for 
these fluctuations can be written in terms of the comoving curvature perturbation 1Z 
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which is gauge-invariant under reparameterizations of time. In practice, and especially 
when carrying the calculation to third order, it is much simpler to work in the ADM 
formalism |32jj, where the metric has the form 

ds 2 = -N 2 dt 2 + hij{ dx i + N i dt)(dx j + N j dt). (21) 

In this representation h V] is the three-dimensional metric on slices of constant t. The 
lapse function N and shift vector N l contain the same information as the metric fields 
$ and B. However, they are chosen in such a way that they appear as Lagrange 
multipliers in the action, so their equations of motion are purely algebraic. After solving 
these constraint equations, the solutions for N and N l can be substituted back into the 
action, thereby avoiding the very lengthy manipulations involved when working with 

GSM2DD- 

All our calculations simplify considerably by working in the comoving gauge, where 
the three-dimensional slices implicit in m are chosen so that the inflaton perturbation 
Sf> vanishes. On slices where Sf> = 0, the three-dimensional metric takes the form 

E3EH§ 

hij = a 2 (t)e 21z 5 ij , (22) 

where the held E has been gauged away by an appropriate choice of the coordinates x\ 
and 7 Z is the non-linear generalization of the comoving curvature perturbation m-\\ 
Although in principle this is a gauge choice, our results will be gauge invariant up 
to reparameterizations of the spatial coordinates. This follows since the quantity 1Z 
is actually gauge-invariant to all orders, being defined by the physical condition that 
comoving observers see vanishing momentum flux {421, 42] • We apply the comoving gauge 
uniformly throughout the present paper. In principle there is some interest attached to 
working with other gauges, such as the spatially flat gauge or uniform density gauge, but 
in such cases the formalism we will describe becomes burdened with a large number of 
extra terms. These terms arise from spatial derivatives associated with inhomogeneities 
which are generically present in <fi, but are absent in the comoving gauge where df = 0. 

3.1. The constraint equations 

With the ADM metric ED, the coupled action 0 reduces to 

A = iy dtd 3 x VhN (i? (3 > + 2 P) + \J dt d 3 x VhN~ l (E i3 E ij - E 2 ), (23) 

where h = det and R ® is the Ricci curvature calculated with hij. The symmetric 

tensor Eij is proportional to the extrinsic curvature of the spatial slices, 

^ij = — N(i\j) (24) 

§ Our notation is chosen to correspond to 32, where the symbol TZ is used for the curvature 
perturbation in the comoving gauge and £ is used for the curvature perturbation in the uniform density 
gauge. In standard inflation these coincide up to choices for signs, but this need not be the case once 
a non-trivial kinetic structure has been introduced into the Lagrangian. 

|| The reader is warned that different conventions for extending (|20ll beyond linear order are employed 
in the literature. The situation is nicely reviewed in ED 
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where | is the covariant derivative compatible with h %] . The N and IV* constraint 
equations are 

r( 3) + 2 p _ 4 XP X - -^r( EijE* - E 2 ) = 0, 


N 2 


N 


{El - ESj) 


- 1 1 3 


= 0 , 


(25) 

(26) 


respectively. 

In solving these equations, we follow and split the shift vector IVj into 
irrotational and incompressible parts, IVj = + IVj, where N ul = 0. After setting 

N = 1 + a, the quantities a, if and IVj admit expansions into powers of 7 Z, 


OL — Ot\ + OL2 T ' ' ' , 


0 = -01 + H-, 

Ni = JV?> + Jvf 1 + - •, 


(27) 


where (for example) a n = 0{lZ n ). We then set the constraints to zero order-by-order. 
The background equation is the Friedmann equation At hrst-order, one hnds from 
the N l constraint that 


«i = dPNV = 0, 


so with an appropriate choice of boundary conditions one can justifiably set 
It follows from the N constraint that^f 


(28) 
= 0. 


= -§ + a- 2 h, ( 29 ) 

where S was defined in (113. As emphasized in E7i. when calculating the action to 
order n in 7 Z, we do not need to compute the order- TV 1 term in N or IV*, since this must 
be multiplying dL/dN or dL/dN' 1 and these are both zero by virtue of the constraint 
equations. In general, one would need all terms up to and including 0(7?. n_1 ), but in 
the present case, terms of order 1Z 2 drop out of the third-order action, so (EHD-dUD are 
sufficient to calculate the order- P? term. 


3.2. The free field action 


Using IHUUD to solve for N and IV* in the action and keeping terms up to quadratic 
order in 1Z, the second-order action is 


S 2 


dr d 3 a: a 2 




e{dlZf 


(30) 


in agreement with the action for ^-inflation calculated by Garriga & Mukhanov |20j . 
where r denotes conformal time, defined by d t — a dr, and a prime ' denotes a derivative 
with respect to r. Conformal time during inflation is given to leading order in slow-roll 
by r = -( aH ) _1 . 


The operator d 2 is the solution operator for the Laplacian, defined by d 2 (d 2 4>) = f. 
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In practice, it is convenient to introduce a rescaled field v = zlZ , where z is defined 


by 


z 2 = 


2a 2 S 2 a 2 e 


H 2 


and the speed of sound c s was defined in 
where the operator A satisfies 

8 2 


A = 


+ c 2 d 2 + — 
Sr 2 ‘z ' 


(31) 

In terms of v the action becomes \ f vAv, 


(32) 


The u-propagator between time To and time r is G v (r, To) = iA _ 1 (T, r 0 ), where to reduce 
clutter in the notation we have suppressed the spatial dependence in A. When expressed 
in terms of Fourier modes, this simply means that 

G ” (T ’ To) = /(!y 

where 



G v (k, t )e _ik ' ( - x_y \ 

(33) 

1 

1 

(34) 


Eq. (1341) is known ns the IViukhnnov eejuntion nnd is ecjuivctlent to /\v — o [SZlEHl 110!- 
In general, this equation for G v is not easy to solve. The effective mass z"fz can 
be expressed in the form 

' = (35 ) 
z t a 

where v is a combination of terms that are linear and quadratic in the slow-roll 
parameters m EH- To obtain an approximate solution in standard inflation, that 
is valid to first-order in slow-roll, terms in v that are first-order in £ are treated as 
constants and terms of 0(e 2 ) are dropped. This procedure only makes sense if the 
time derivatives of the 0(e) quantities may be neglected along with the 0(e 2 ) terms, 
which requires |e|, \rj\ <C 1. Of course, this is nothing more than the familiar slow-roll 
approximation of standard inflation. However, in the present context there is an extra 
condition arising from the requirement that c s must also be kept constant" 1 ". The error 
arising from this latter approximation will be at least as significant as that arising from 
the mass term, so consistency requires that s = 0(e 2 ). After taking into account Eq. 
(Hi, this implies that u must itself be of order e. In other words, the approximate 
solutions of (El are only valid if c s is sufficiently close to unity, to within a quantity 
that is small in the slow-roll limit. 

To first-order in this generalized sense of slow-roll, the 7Apropagator satisfies 

(1 — i/cc s r 0 )(l + i kc s T)e~ lkCs( ' T ~ T °' > r > r 0 
(1 + i/cc s 7o)(l — i/cc s r)e lfcCs( ' r_ro ' ) r < To 


n , s H 2 1 
G R (t,T 0 ) = —- X 


(36) 


+ The much more complicated case where c s may have some appreciable evolution was considered in 

E31I221 
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where we have chosen boundary conditions so that Gn behaves like the flat space 
propagator at very early times, when the mode is deep within the horizon and cannot 
feel the curvature of spacetime. This corresponds to the Bunch-Davies vacuum H- 
The power spectrum of TZ is easily obtained from © and was derived in [20] • On 
large scales, dropping less singular pieces as k —> 0, the two-point function is* 


(7e(k 1 )7e(k 2 )) = (2vr) 3 ( 5(k 1 + k 2 )P(k 1 ), 


Tt2 1 

p < k '> = 4FII’ 


(37) 


where k — |k|. To turn this into a power spectrum one takes the coincidence limit to 
find the dispersion, a 2 = (TZ(x) 2 ), and then evaluates its logarithmic derivative: 

da 2 1 H 2 


A ^ d hi k 8tt 2 £ 
The tilt of this spectrum is given by 

din A 2 (A;) 


n 


1 = 


din k 


2e — q. 


(38) 


(39) 


The quantity TZ is conserved outside the horizon HUES!, which is the analogous 
result to that of conventional slow-roll inflation mm- 


4. The third-order action 


4-1. General form of the action 


We now turn to the central calculation of the present paper, a determination of the third 
order piece in the coupled action (TTtI) . In order to compute this, we need expressions for 


R^ 3 \ P(X ) and — E 2 . It is easy to show that R ® = —2a 2 e 21l [(dTZ) 2 + 2d 2 TZ], 


To calculate P(X, 0), we use the fact that 

i iz iz 2 IZ TZ 3 

—- = 1 — 2-2o 2 + 3— + 6a 2 -4—. 

N 2 H 2 H 2 2 H H 3 

Since we have chosen a gauge in which <50 = 0, this means 


(40) 



„k 2 , n 

2a 2 + 3-yr + 6a 2 — 
ri 


Ai 

H 5 


P 


x 


PZ 2 TZ TZ 3 

4— + 8a 2 -12— 

H 2 H H 3 


Pxx ~ ~^X 3 P 


3!' 


,xxx 


TV_ 

w 


(41) 


No derivatives of P with respect to 0 occur, since the inflaton field takes its unperturbed 
value. This is one convenience of working in the gauge 5<f> = 0. Finally, using the 


* Our expression for the power spectrum differs by an overall factor of c s 1 from the corresponding 
expression in EDI- At leading order in slow-roll this is harmless, because we have already seen that 
c s must be equal to unity to within O(e). However, if one wishes to expand consistently in the small 
parameters of the approximation, one should set c s = 1 exactly to leading order in l.‘17li - (See |2.'1| for 
a more detailed discussion of this point.) 
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connection derived from h %3 , and our choices for N and N t , we obtain 


E ij E tj — E 2 = -6 (H + 1l) 2 + 


AH 


iz. 


(1 + ^)e- 2n (d 2 fj + dlZ ■ df> + N ■ dlZ) 
H 


-^-e- A7Z (d' 2 iljd 2 f + 2 d 2 fj[dlZ ■ dip]) 


+-re 
a 4 


—47?. 


Vij Vj + 2 ipjjNij - 4 ‘Rj'ip'ijipj + 2d 2 ip[&TZ ■ dip]). 


(42) 


Once one has collected terms and integrated by parts where possible, it turns out 
that all terms involving ip 2 and Ay ; either cancel among themselves or reduce to total 
derivatives, which can be discarded. The remaining second-order contribution from the 
ADM quantities N or N’ 1 is just proportional to a 2 , and can be written as 


S 


ADIVr 


1 r drd 3 x AHald 2 ^--^d 2 Tl + 


;E1Z 


(43) 


2 J ~ ri H~ ' H 2 

This vanishes when ipi takes its on-shell value El, since it is proportional to a 
constraint. (This justifies the statement made in Section ETT1 that it is only necessary to 
calculate N and N t to first-order.) It is most economical to rewrite the resulting action 
in terms of the quantities £ and A of (hid- After integrating by parts, discarding total 
derivatives, and using the background equations of motion, the third-order contribution 
to the action can be written in the form 


S 3 = - J dr d 3 x a 3 


2 H 


n 3 


jin 2 


i 2 Hi nmf-( 2 E + 4 A)^ 5 + 6 E : ff2 


-VVi/v.vi., 

a 4 


1 n 


1 n 


- —J 19 fa 9 ^ + —pH + 

a a ri a a n 


(44) 


To complete the reduction, we must replace ipi with its value, given by Eq. HI). 
In doing so, it is very convenient to make use of the equation of motion derived from 
the free, Gaussian theory defined by (EDD- If we introduce a new quantity A, satisfying 
^2 


a = 

the field equation which follows from m can be written as 
5L 


511 


dA ZJ A 

= Tt +HK 


sd TZ. 


(45) 


(46) 


This vanishes when 1Z is a field mode which solves the equation of motion (El of the 
Gaussian theory, but 5L/51Z |i will be non-zero when 1Z satisfies the equation of motion of 
the full interacting theory that takes into account the 1Z 3 vertex that we are computing. 
By using (TIBI) , we can trade derivatives of A, which will appear after integrating by 
parts in dH , for simpler terms involving HA, 0 2 1Z and 5L/51Z\\. Eventually, we will 
be able to remove the 5L/51Z\i terms by a field redefinition, resulting in a considerable 
simplification of the 1Z? vertex. 

Proceeding in this way, one finds that the action can be written in the form 


S 3 = -J drd 3 x a 3 1 


2 iz 3 

- —£lZ(dlZ) 2 - (2£ + 4A) — + 6 E1Z1Z 
a z 
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44 (d'R-) 2 ~ 4 An ,iX,i + 4 Aen ,iX,i + 


a 4 H ' a 
+ terms involving 5L/A7?.|i^ , 

where the terms involving SL/6R\i are given by 


2 a 


(47) 


a'H- 





+ (d 2 ^X,i + A77 j) ft? ' 




(48) 


and we have set x — d~ 2 A. These terms will not be important in what follows, because 
they all contain at least one derivative of R and therefore vanish outside the horizon, 
where R approaches a constant. When determining the order in slow roll of each of 
the terms in m-m, it should be noted that y = 0(e): x (and A) are first-order in 

slow-roll. 

After further integrations by parts and use of the Gaussian held equation © , we 
find that the entire three-point vertex can be expressed in terms of the how parameters 
e, 77 , Ex, u and s: 


S 3 = - J df d 3 a; a 3 

2 e 


a 2 cl 


4 / 1 e 

-£ --- S + U 

3 l 3c 2 Ex 


(e — 2s — uc 2 s )1Z(d1Zy 



1Z1Z 1 


44 ^ n ,iX,i 


a 2 cl 


K 2 k £ 4 (1 
cl d t\ci 


£“ _ L 9 1 


— K-Tl + 


(49) 


Although this has the appearance of a series expansion in terms of how parameters, this 
expression is in fact exact to 0(7Z 3 ), given that interactions of both gravity and the 
scalar held with the other contents of the universe have been neglected. This vertex 
should be supplemented by terms proportional to 5L/S7Z\\: 


S- 


3, Gaussian 


1 

2 


' dtd z x 7 l^n 2 


8L_ 

Jiz 


+ 


(50) 


where we have omitted the derivative terms in (HHD, which vanish outside the horizon. 

Each of the parameters appearing in (BHI) have comparatively simple interpretations. 
The parameters e and r] describe the non-gaussianity produced by coupling a scalar held 
(with any given kinetic structure or self-interaction |27j ) to gravity. Terms containing 
u measure how far the sound speed deviates from the speed of light: in other words, a 
dispersion relation different from the canonical case E 2 = p 2 + m 2 acts as a source for 
non-gaussianity. There is another source of non-gaussianity which arises from changes 
in the speed of sound, as measured by s. This particular source appears in combination 
with e/ex, which measures the non-trivial nature of the kinetic structure in P(X, <f>). 
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We only calculate to leading order in slow-roll. From a practical point of view, this is 
necessary given the complexity of the calculation beyond leading order, but since the 
predicted non-gaussianity is rather small E3. only the leading-order effect is ever likely 
to be observed. Restricting (HID to leading order terms - in this case, they are 0 (£ 2 ) - 
we find that the three-point vertex is well-approximated in the slow-roll regime by the 
expression 


S 3 = 


d t d 3 x a 3 


4 

~ 3 £ 


u + 


£ s\ 1Z 


£x3/ H 


+ 2e{3 u + e)V,il 2 


+-^£(£ - u)TZ(dTZ) 2 - -4 eilTZ } iX,i , (51) 

cr cr J 

together with one supplementary term proportional to 5L/87Z\i, as given by Eq. (15011 . 
This is also a leading-order term in slow roll. Despite appearances, the term involving s 
really is of order £ 2 , since we have already seen from Eq. (USD that £x = 0 (£ 2 ) whenever 
a slow-roll regime exists. Although in principle we could proceed to the calculation with 
this vertex in its present form, it is worth performing a further integration by parts to 
remove the term involving x,i- After doing this, and using the Gaussian held equation 
(USD , we can rewrite (ED as 


S 3 = / dtd 3 x 


£ S 


n 3 


- -a 3 £ (u +-— + 4a 5 £(2u + e)HTZ 2 d^ z TZ 

£x 3 / H 


4 a 6 £nd z nd- 2 n 


(52) 


together with some terms that are proportional to the Gaussian equations of motion, 

(53) 

where 


S: 


3, Gaussian 


df d 3 x aT 


5L 

PR 


t = (77 - u - £)n 2 + 2£d~ 2 {nd 2 n) ( 54 ) 

includes the contribution of (ED- There are also terms that are of order 0(£ 3 ) and 
higher, as well as terms which are proportional to the equations of motion, but these 
are irrelevant since they contain derivatives of 7 Z. 

Finally, we should confirm that the form of this vertex reproduces known results 
in standard, single held inflation. In fact, assuming the canonical form P(X, <f>) = 
|A" — V(4>) for the polynomial P, it is easy to show that (I521) - (|5H) . reproduce the 
corresponding results of EH- 


4 . 3. Field redefinitions 

Before describing the calculation in detail, we hrst show that a held redefinition of the 
form 1Z I—► lZ n + F(lZ n ) can be used to eliminate terms proportional to the Gaussian 
equation of motion, 5L/8TZ\i. In principle we could evaluate these terms with the others, 
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but the most economical way of accounting for them is to remove all such terms via a 
suitable field redefinition, and then incorporate their effect into the correlation function 
by applying Wick’s theorem. Any terms which vanish outside the horizon, such as 
those omitted in © or (H3D, will make no contribution to the correlation functions on 
superhorizon scales. Hence, they may be legitimately ignored. 

A field redefinition of the form TZ > 7 Z n + F(JZ n ), where F is quadratic in lZ n , 
has no effect on any of the 0(7Z 3 ) terms in On the other hand, its effect on the 

Gaussian term m is to transform 


S 2 [K] 



^e[nl + 2n n F(n n )} 


After integrating by parts, this is equivalent to 


a e{dKl + 2dK n dF) 


(55) 


(56) 

According to this general argument, we need to make a held redefinition 7 Z i— > TZ n + \Fhi 
order to remove the terms arising in Eq. (1531) . Once this redefinition has been made, we 
can calculate the three-point function corresponding to the vertex 0-0 by working 
only with Eq. (E2D rewritten in terms of 7 Z n . 

In the following section, we employ m to calculate the three-point function 
(7^(k 1 )7^(k 2 )72.(k 3 )) of the fluctuation. 


S 2 [TZ] i-> j S 2 [JZ n \ -2 J dtd 3 x Fa — 


5. Path integral formalism for the three-point function 

The three-point function we wish to calculate is (7£(ki)7£(k2)7£(k 3 )), which measures 
correlations produced by the vertex H52D-0- In this section we assume that the held 
has been redefined to remove terms proportional to the equation of motion, but for 
clarity in formulas we drop the subscript l n\ This expectation value is to be taken in 
the interacting vacuum |G) of the theory, that is, 

(7^(k 1 )7^(k 2 )7?.(k 3 )) = (fi|7e(k 1 )7e(k 2 )7e(k 3 )|fi). (57) 

5.1. The interacting vacuum 

In order to assist the comparison of our formulae with those of nam], we begin by 
briefly reviewing the standard construction of the interacting vacuum, following B3- 
The quantization of a free held such as m, corresponding to the Gaussian action m, 
proceeds by writing 

TZ(t, k) = 7Z C \(t, k)a(k) + TZ* cl (t, k)a^(k), (58) 

where a(k), a^(k) are annihilation and creation operators in the usual fashion, and 1Z C \ 
and 7Z* cX represent a classical solution of the held equation, such as that given in Eq. 
and its complex conjugate. The free vacuum |0) is constructed so that a(k)|0) = 0 for all 
k. After introducing a self-interaction, such as the 1Z? vertex (I49[l . the mode functions 
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of 1Z can no longer be calculated exactly because the interaction mixes Fourier modes. 
This implies that a construction such as that given in Eq. © is no longer possible. 

As a result, we must construct the interacting vacuum in a different manner. One 
approach is to begin with (USD when the theory is approximately non-interacting, and 
evolve it forward according to the familiar Heisenberg rule 

K(t, k) = e iH{t - to) n(t 0l k)e~ iH{t - t0 \ (59) 

where t 0 is a fixed fiducial time at which (EHD was constructed, and H is the Hamiltonian. 
For example, the principle of cluster decomposition usually means that we can construct 
© as a solution of the Gaussian theory at asymptotic past infinity. In the cosmological 
context, this is when the mode corresponding to k is deep inside the horizon. 

The Hamiltonian H can be split into a piece Hq, corresponding to the Gaussian 
action ©, and a piece Hj, corresponding to the self-interaction (14911 . It is then 
straightforward to verify that the quantum operator TZ satisfies 

K(t, k) = U\t,t 0 )Mt,k)U(t,t 0 ), (60) 

where the interaction-picture field 7 Zj is a solution of the free held theory and the time- 
evolution operator U is defined by 

H = Texp(-i£dCtf/(C)), ( 61 ) 

where T represents the time-ordering symbol. The interacting vacuum Q) should be 
destroyed by annihilation operators corresponding to the full interacting theory, not the 
operators a(k) of the Gaussian theory. To obtain |fi), one evolves |0) for some time T, 
such that 


e- iHT |0> = e" iEoT |fi)(fI|0) + ^ e~ iE " T \n)(n\0), (62) 

nyo 


where Eq = (0|i7|0), and the {E n } are the spectrum of Ed. Since |fi) is the vacuum of 
the theory, it follows that E 0 < E n for any n, so a slight rotation of T into an imaginary 
direction, T —> oo(l — i<5), implies that all terms from the sum over n / 0 become 
exponentially small when compared to the term involving |Q). It follows that Q) can 
be written 


P) 


lim 


,-LffT 


T-oo(l-i<5) e~ iE ° T (fl |0) 


10 } ■ 


(63) 


Combining this expression with the expression for the Heisenberg held, Eq. m, implies 
that one can compute the correlation functions of the interacting theory according to 
the rule 


(n\Tn(x) ■ ■ - n( y )\n) 


lim 

T— xx)(l—i<5) 


(0\TR.!(x) • • ■ My) exp (-i /5 T d t Hj(t)) |0) 
(01 exp i jl T dt Eff(t)') |0) 


(64) 
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Depending on the details of how |Q) relates to |0), one may wish to choose a more 
general contour of integration in m to project out the true interacting vacuum. Using 
this prescription, we can calculate the three-point function of the interacting 77 3 theory 
from the usual path integral rule: 

/ [d77] 77(f ) 3 exp(i f c L) 


(D|77(f) 3 |D) = 


(65) 


/[d77] exp(i/ c L) 

where C is a contour of integration, chosen to select the interacting vacuum as described 
above, and S = f c L is the action. For the present theory this consists of the Gaussian 
piece HD, which we write as S 2 , and the three-field vertex which we write as 

S 3 . In principle, we should include diagrams containing an arbitrary number of loops 
when evaluating (ESI), but because S 3 is one order higher in slow-roll than S 2 , the loop 
expansion is effectively the same as the expansion in slow-roll. As a result, if we are 
only interested in the leading-order slow-roll dependence, we need only calculate to tree 
level. 

After expanding the interaction part of e lS as a series in the usual manner, and 
recalling that the integration of any odd number of 77’s against a Gaussian measure is 
identically zero, we find that the tree-level three-point function is given by 


(D|77(f) 3 |D) tree = 


J[d77] 77(f) 3 (/ c L 3 )exp(i/ c L 2 ) 


( 66 ) 


/[d77] exp(i f c L 2 ) 

As an explicit example, let us consider a vertex of the form L 3 = p77 2 cU 2 77 
(where g is some coupling constant) which arises in (15 1 1) after going to conformal 
time, and which we will need in the following section. In this case, one finds, after 
performing the functional integrations, that the three-point function can be written as 
(77(ri, y 1 )77(r2, y 2 )77(T 3 , y 3 )) = i g f c dr]d 3 x Q, where the integrand Q satisfies 

G = J2 j-( 7 ^( r i,yi)77(^,x))-^(77(r 2 ,y 1 )77(p,x))c)" 2 -^-(77(r 3 ,y3)77(p,x)) (67) 

pairings dr l dr l d V 

and the sum is over all ways of pairing 77 (ti, y 1 ), 77 (t 2 , y 2 ) and 77 (t 3 , y 3 ) with fields in 
L 3 . Notice that in principle, the path integral sums over all pairings of the 77’s among 
themselves, including pairings of external fields with fields at the vertex, and pairings of 
external fields among themselves and vertex holds among themselves. However, there is 
no contribution from pairings where two external fields or two internal fields are paired, 
because the resulting amplitude is proportional to the tadpole for an 77 to emerge from 
the vacuum. Since we are assuming that the vacuum is stable, all such amplitudes 
vanish. This prescription is equivalent to dealing with only the connected three-point 
function. 

After going over to Fourier space and using Eia) for the propagator, one finds that 

/ H 2 \ 3 /(1 + ntesT ) 3 


(77(k 1 )77(k 2 )77(k 3 )) = -u/(27r) 3 <$(5> 

i 

(1 - ikc s T ) 3 


. 4 £C S 


n ,kf 


2 /„ 2, n 2 o~ lKc dV— T ) 


dr} kfkpq 


n ikf 


+ perms, 


dp kfk^r] 


2 h 2 n 2 P iKc s(.v-P 


( 68 ) 
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where we have evaluated the correlation function at equal times Ti = r 2 = r 3 = r and 
K = k\ + fv 2 + &3 is the total momentum. The mode is deep inside the horizon in the 
far past and oscillates rapidly. In this regime, there is no contribution to the integral 
once we have rotated to Euclidean space. In the far future, 1Z tends to a constant and 
there is also no contribution. The dominant behaviour is determined by the modes’ 
characteristics around horizon crossing. Therefore, even though the propagator CSD 
is only valid under the assumption that u, £ and g are small, Eq. flUED will still be a 
good estimate if the slow-roll parameters are sufficiently small around horizon crossing 
and H , e and g are evaluated when these /c-modes crossed the horizon. Taking these 
considerations into account, we arrive at the simple formula 

{K(k 1 )K(k 2 )TC(k 3 > = -i 9 (2x) 3 i(X:k i )iJll ^ f^inhlkW^ + c.c., (69) 

plus permutations, where ‘c.c.’ denotes the complex conjugate of the preceding term, 
and we have taken the asymptotic value by sending r —> 0. 


6. The scalar non-gaussianity 


6.1. The three-point function 

The interaction vertex dS2D contains three terms. Two of these, proportional to 1Z 3 
and 7Zd 2 TZd~ 2 7Z, are new and arise because the effective speed of sound c s ^ 1. They 
vanish in the limit where c s = 1. The third, proportional to the TZ 2 d~ 2 7Z interaction, 
generalizes the corresponding term derived in [2Zj to the case where « ^ 0. At the 
tree-level approximation, the contribution that each of these interactions makes to the 
three-point correlation function is expressed in Eq. m- We now proceed to evaluate 
each of these contributions in turn. 


1Z 3 interaction. Using the propagator (1361) . and evaluating this vertex in a similar 
way to that outlined above, we find that this term gives a contribution 
H 4 u-\- es/?>£x 


i(27r) 3 5(^k, ; 


klk^kfrf dr/ e lKcsV + perms + c.c. 


(70) 


. 3 - 2% 2 Ylik? 

Since there are no poles, there is no obstacle to rotating the contour of integration 
so that it lies along the imaginary axis. This can only be done in one direction, 
the direction compatible with the infinitesimal displacement T —> oo(l — id), as 
described above. This means that we must rotate C clockwise onto (oo, 0)i, thereby 
yielding a total contribution equal to 

H 4 u + £s/2>£x k\k\k\ 


- (2vr) 3 d(^k f 


3 • 2% 2 n, kf K 3 


+ perms + c.c. 


(71) 


lZ 2 d 2 TZ interaction. The second type of interaction contained in the vertex is the 
u / 0 generalization of the equivalent term that arises in the case of canonical 
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inflation. This is the vertex we calculated in detail in (I(i9l) . with g = da 5 £(2u + E)H. 
The contribution is therefore given by 

- i(27r) 3 h(^k i )^— — \j(2u + e) [ a 3 r) 3 kfkl dr? e lKcsV + perms + c.c. (72) 

i 2 £ Ui K JC 

Since ag = —H^ 1 during inflation (to leading order in slow-roll), all g terms drop 
out of the integral except those arising in the exponential. After rotating the 77 - 
integral to the imaginary axis as in the previous interaction, we find that this 
contribution is given by 

(2tt) 3 <5(^ k,) ^ 2 ^3 ( 2m + e) 1 ^- + perms + c.c. (73) 


TZd 2 lZd 2 1Z interaction. The remaining new piece in the vertex involves the term 
Tid 2 7Zd~ 2 1Z. This yields a correlation of the form 
H 6 


-i( 27r ) 3 5(^k ?: )- 


u 


j a 2 r] 2 klkl(l — i k 2 c s 7j)e lKcsT1 + perms + c.c., 


(74) 


2 % 2 n M 

with our usual conventions for K. After rotation of the contour C, we hnd that this 


contribution reduces to 
H 4 


- ( 2 T) 3 <KH k * 


u 


2 4 e 2 [I; kf V K K ‘ 2 


k 2 k 2 k 2 k 3 N 

rv 1 run Iv 1 A/O 

+ -rrf- + perms + c.c. 


(75) 


It is now necessary to take into account the contributions to the correlation function 
which enter via the field redefinition that was introduced to remove the terms in Eq. 
(j5Up. There are effectively two redefinitions, as summarized in Eq. (ED, and we consider 
each of these in turn. 


• The redefinition 1Z 7 Z n + (7 j — u — e)lZ 2 a / / f. Under a field redefinition of the form 

1Z 1 — 7 7 Z n + qlZn, Wick’s theorem guarantees that 

{7Z(x)7Z(y)7Z(z)) = q(7Z(x)7Z(y))(7Z(x)7Z(z)) + perms. (76) 


Any correlation function of the form (lZ(x)TZ(y))(IZ(x)TZ(z)) can be written in 
Fourier form by taking appropriate transforms in x, y and z. After making the 
necessary transformations and integrating out the (5-functions (CUD, we obtain the 
product of two copies of the spectrum, 


-(77 
4 V 


u 


H 4 

16 E 2 


d 3 h d 3 k 2 


pix (ki+k 2 )g—iy kig—iz k 2 _ 


+ perms, 


(77) 


(27t) 6 kfk 2 

I 11 order to express this result in a more familiar form, we introduce an auxiliary 
integration to represent the combination ki + k 2 , constrained by a (5-function, such 
that the contribution is equivalent to 

1/ N H 4 , 

4^-«- £ )^( 27 0 


d 3 ki d 3 k 2 d 3 k 3 

( 2 Iff 


<5(^ kOe ^ 3 '^ 1 y+k2 z) - 3-3 + perms. (78) 
i k 1 k 2 


In this form, the fc-integrals and appropriate factors of 2n can be stripped off to 
yield the correct /c-space contribution to the correlation function. 
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The redefinition 7Z i—> lZ n + (e/2)d 2 {JZ n d 2 lZ n ). In this case we obtain contributions 
to the three-point correlation function of the form 

£ 

2 * 

Introducing the Fourier transform of the two-point functions, as outlined above, 
and integrating out h-functions, it follows that this term must be identical to 


(K(x)H(y)K(z)) = W 2 ((K(z)Jl( y )%(K(x)K(z))) . 


(79) 


H 4 

—, d~ 2 
32e 1 . 


d 3 ki d 3 k 2 


c ix (ki+k 2 ) e iy ki e iz k 2 _ 


+ perms. 


(80) 


(27t) 6 Ayfcf 

As before, by introducing an auxiliary constrained integration to represent ky + k 2 , 
this can be written in the form 

H 4 1 




2 • 2 4 e k\ kfkf 


+ perms, 


(81) 


after /c-space integrations and factors of 2 n have been removed. 

Bringing all of these separate contributions together, we arrive at an expression for 
the three-point scalar correlation function: 


(7e(k 1 )7e(k 2 )7e(k 3 )) 


(2vr) 3 5(Ek i ; 


H 4 


2% 2 n, kf 


A, 


where the ^-dependence is determined by the function A: 

A = ^(u + e) E k- k 2 ~^{u + —0 k\k\k\ 

i>y 


2 u 
K 2 





(82) 


(83) 


In evaluating this expression, we have assumed that Ay ~ k 2 ~ k^, so that the epoch of 
horizon crossing for each of the fc-modes is comparable. Eqs. constitute the 

principal result of this paper. 


6.2. Consistency conditions 

The full three-point function 1S2NH3J is rather complicated, and one would like to 
have an independent check of its validity. This may be achieved by relating the three- 
point function to the spectral index of the two-point function in an appropriate limit, as 
discussed in |2Zl EH US!- We begin by briefly reviewing the argument of P within the 
context of models where c s 7 ^ 1. The aim is to write down a long-wavelength gravitational 
consistency condition between the two- and three-point correlation functions. In the 
limit where one momentum, say k 3 , is much smaller in magnitude than the other two, 
such that C fci ~ b, the remaining momenta ki and k 2 become approximately equal 
and opposite. From the point of view of physical perturbations, this means that the 
wavelength of the k 3 mode has been made arbitrarily long. A mode such as this, which is 
deep in the infrared and far outside any individual observer’s horizon, is effectively just 
a renormalization of the background theory, which in this case corresponds to a rescaling 
of the spatial coordinates x l by a factor e ^ 3 ~ (1 + 7 Z 3 ), i.e., 5x = IZ 3 X, where 7 £ 3 is the 
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amplitude of the k 3 mode. Such a rescaling implies that the wavenumber of any mode 
still inside the horizon must be compressed by a corresponding amount 5k = —77 3 /c and 
this compression results in a change in the epoch of horizon exit, kc s = aH, such that 


5k ■ c s — 5a ■ H. (84) 

If we use the standard relation 5a = aH5t, this implies that 5t = —1Z 3 /H 
independently of the value of c s . As a result, we expect the three-point function in this 
limit to express how the k 3 mode correlates with the change of the two-point function 
of ki and k 2 induced by differences in the epochs of horizon exit, i.e., 


(7e(k 1 )7e(k 2 )7e(k 3 )) - (27 t) 3 5(^kj)(l — n)P(ki)P(k 3 ), (85) 

l 


where P(k) was defined in (liTTI) . and n is the spectral tilt (|39l) . This consistency equation 
was rederived on kinematical grounds in inn, where it was emphasized that <|B5D holds 
independently of the details of the matter Lagrangian. Therefore, one should expect 
dH2D-(jHSD to obey dHSD as | k 3 1 is sent to zero. 

In fact, it is straightforward to verify that the momentum sums in (1851) behave 
according to 


h 2 h 2 
\ - G K'j 

if.-jr 

i>3 


v 2) 


E 


K 2 



E^ 3 -2 kl 


E k i k j ~ 2k h 


( 86 ) 


whereas kfkfkf 


0. This behaviour is sufficient to show that 

H 4 1 


(n(k 1 n(k 2 )n(k 3 )) - ^(E^ 


2 4 e 2 k\k 3 


(2 e + v), 


(87) 


which reduces to m when the constant of proportionality is rewritten in terms of the 
scalar spectral index. 

An estimate of the level of non-gaussianity generated purely by the scalar held, 
rather than its coupling to gravity, can be derived by specifying £ — r) — 0. In this 
regime, the ^-dependence of (1831) reduces to 


k 2 k 2 

A\ £=v =o = 4 u ^2 

i>j 


1,2 1,21,2 
, IX /1 A/Q r\j Q 

_ Ail —i —-—— 

K 3 


_ 7 / _ 

- 2 -E^p-fE^ 

*7 1 


( 88 ) 


and, in particular, A \ £=v=0 —> 0 as k 3 —> 0. The behaviour of A\ £=ri= Q in this limit 
guarantees that UEZD is independent of the speed of sound. 


6.3. Observational considerations 

The level of accuracy currently achievable by observations is insufficient for a projection 
of the three-point function onto the CMB sky to be measured directly pen. 

Instead, observational limits are set on a parameter / NL |4|| 4j2j, which measures the 
departure from Gaussianity of the gravitational potential. This potential IS3, together 
with other detailed physics at the last-scattering surface, controls the character of the 
CMB anisotropies we measure at the present epoch. These effects must be taken into 
consideration when calculating the non-linearity /nl that a given model of inflation will 
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imprint in the CMB inn E2J- Present-day observations impose a conservative upper 
bound of |/ NL | < 100 J53], whereas a non-gaussian signal exceeding / NL > 5 should be 
detectable from forthcoming observations of the Planck satellite |48j . 

Quite generally, one can write /nl as the sum of a superhorizon piece, which is 
is produced by gravitational evolution of the superhorizon perturbations after inflation 
has ended and a primordial piece, fff L , related to the non-linearity in the curvature 
perturbation 77, which should properly be regarded as a kind of initial condition. This 
initial condition is most conveniently characterized by writing 

n(x) = n g (x) - 5 /r * ( n g (xf - (nl)), (89) 


where TZ g {x) is a Gaussian field, and a convolution product * is involved since in general 
may be a function of scale. The standard single-field inflationary scenario predicts 
/nx ~ 0(e) |22, !, whereas the universal superhorizon contribution is ~ 0(1). In 
consequence, a positive detection of primordial non-gaussianities from WMAP or Planck 
would represent one possible way of ruling out such models. We will see that the 
prediction for non-canonical inflation is similar. 

The field (I89j) is just a particular example of a redefinition of the form US) , so it is 
straightforward to express /nl in terms of A: 


J NL — r 


5 A 


(90) 


6 El k? ’ 

Eq. (1HU1) can be written explicitly in terms of momentum sums if desired, but the result 
is complicated and ultimately not informative. A better estimate of the size of any non- 
gaussianity produced in this model is given by evaluating fff L on equilateral triangles, 
where k\ — k? — k-j. It is a simple matter to calculate 


/, 


n 

NL 


equilateral 


275 

972 


u + 


10 e „ 55 , 

729£^ S- 36 £ ~ 

0.28m + 0.02—s - 1.53e 

£x 


12 V 

- 0.4277. 


(91) 


This value of /^ L agrees with the corresponding result of Maldacena |2IJ for 
standard inflation, for which s = u = 0, and can be compared with predictions for similar 
quantities in alternative models pad. For example, Creminelli m has considered a 
model that in our notation corresponds to specifying 


P(X, <t>) = X/2 - V(4>) + x 2 /sm\ 


(92) 


where M is some mass scale. As discussed in the introduction, M should be associated 
with a large renormalization scale which cuts off the details of the physics in the far 
ultraviolet limit. Although gravity was neglected in the analysis of d, the scale M is 
related to the speed of sound in such a way that 
2 M 4 + X 


(? = 


2 M 4 + 3 A 


( 93 ) 







Primordial non-gaussianities in single field inflation 


23 


Moreover, to leading order in M 4 , the slow-roll parameters u and s can be written as 

2X X 

U =- r*j - 

2 M 4 + X M 4 

e 12 XM A X , . 

£ X (2M 4 + X)(2M 4 + 3X) M 4 ' { } 

Thus, u is small and slow-roll (in the generalized sense of Section o> is a good 
approximation if M 4 X, or equivalently, if \f\ -C M 2 . In this limit, and setting 
e — r] — 0 to decouple the contribution from gravity, we obtain a non-gaussianity of 
magnitude 


fK _ 

J NL — 


35 4> 2 
108 M 4 



(95) 


in exact agreement with M- Although the non-gaussianities can be made arbitrarily 
large by sending M 4 to zero, it is important to note that u ceases to become small in 
this limit, and indeed c s —> 1 /-s/3. Therefore, the slow-roll approximation will eventually 
break down well before the non-gaussianities have become large. 


7. Conclusions 

In this paper, we have derived the three-point correlation function of scalar fluctuations 
generated during an early inflationary epoch in the case where the inflaton Lagrangian 
is an arbitrary function of the held and its first derivative. This has applications to the 
^-inflationary scenario USED] and inflation driven by a ghost condensate 00. More 
generally, it applies to any theory containing higher derivative operators which descends 
from a higher-dimensional supergravity or superstring compactihcation, or includes the 
effect of radiative corrections. 

Our computation reproduces previously derived results nam in the appropriate 
limit of a canonically coupled scalar held, and respects the consistency relation proposed 
on kinematical grounds by Maldacena J2ZJ and Creminelli & Zaldarriaga [31 J- In 
the latter case, the reduction of the three-point function (1%21) - (1H31 ) to the consistency 
condition ( 0 ) is non-trivial, since the factors of u which are distributed over the various 
/c-dependent terms must sum to zero as k 3 —> 0. 

This pattern of fc-dependence in the three-point function is quite different to 
that of a slowly rolling held with a canonical kinetic term, but similar to results 
derived in previous attempts to go beyond canonical inflation 0 0. In principle, 
this modification in angular dependence provides a sharp discriminant between these 
different classes of models. However, the overall scale of the non-gaussianity m is small, 
being proportional to a sum of how parameters which are suppressed below unity when 
slow-roll is valid, as is required for self-consistency of the calculation. This implies that 
if the density huctuations which seeded the CMB anisotropies were generated during an 
epoch of slow-roll inflation of the type we have considered, the level of non-gaussianity 
will be small even in the presence of higher-derivative operators of the inhaton held. 
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If the slow-roll conditions described in Section o do not apply, then the situation 
is less clear. It is possible to have quasi-exponential inflation, e -C 1, even if u 
or s are large. In this case, a relatively large measure of non-gaussianity could 
be produced during inflation by reducing the speed of sound, through a suitable 
choice for the functional form of P(A", </>). This possibility may become attractive if 
forthcoming observations indicate the presence of an appreciable measure of primordial 
non-gaussianity in the CMB. However, whenever u is large (so c s departs from the speed 
of light by a significant amount), the slow-roll approximation begins to break down and 
the calculation we have described is no longer appropriate. In this limit, we expect the 
reduced speed of sound to cause the spectral index to deviate from unity. Consequently, 
a larger non-gaussian signal will typically be accompanied by a larger deviation from a 
Harrison-Zerdovich spectrum. 

To make further progress, it will be necessary to obtain the solution of the 
Mukhanov equation when one or more of the slow-roll parameters undergoes significant 
variation. Some progress has recently been made in this direction 1221123, but in any 
event, the effective held theory for the inhaton will probably cease to be valid before u 
can be made very large [Ill- 
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